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Abstract. We construct stable configurations of n overlapping discs 
of radius r in a unit square, with r = 0{l/n). By a result of Diaconis, 
Lebeau, and Michel, this result is best possible, up to a constant factor. 
A consequence is that the Metropolis algorithm, a well-studied Markov 
chain on the hardcore model, is not irreducible in this range of parame- 
ters. 



1. Introduction 

Sphere packing problems have been studied for at least the past few cen- 
turies, since Kepler and Newton. The densest packing of spheres in three 
dimensions is easy to guess and must have been known for quite a while, 
and was conjectured by Kepler, but the fact was only proved recently by 
Thomas Hales [[lo|] . (Under the assumption that the centers of the spheres 
lie on a lattice, the problem is much easier, and the answer was already 
known to Gauss.) Hales' proof famously relies on a difficult computer cal- 
culation. 

Even the two-dimensional case of sphere packing, the densest packing of 
circles in the plane, is nontrivial. This problem was solved by L. Fejes Toth 
only in 1940. Since then, there has also been a wide literature on densest 
packing of spheres in a bounded convex region, such as a circle, triangle, 
or square 12; 12]. The case of a region with a boundary seems much 



harder in general. 

We use the words "disc packing" in this article to distinguish from the 
area of circle packings and discrete conformal geometry [15]. As an inter- 
esting historical note on disc packings, in 1964 Boroczky used a subtle con- 
struction to disprove a conjecture of Fejes Toth, that a stable arrangement 
of discs in the plane, one where each disc is held in place by its neighbors, 
must have positive density [j^. Boroczky 's counterexample is at the heart 
of what we do here. 

Regular 12-gons and equilateral triangles are well known to tile the plane 
vertex-transitively [9]]. If suitably scaled discs are placed at the vertices 
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Figure 1 . A tiling by regular dodecagons and equilateral 
triangles, and placing suitably scaled discs at all of the ver- 
tices, results in a stable configuration in the plane. 

of the polygons, one obtains a stable configuration, as in Figure [TJ Ap- 
parently some have conjectured that this configuration has minimal den- 
sity [,16,1 . One computes that these discs asymptotically cover a fraction 
of (7\/3 — 12)n ^ 0.390675 of the plane. We note another planar stable 
configuration in Figure[2l which we discovered while preparing this article. 
This planar configuration has a slightly greater density of 0.4569, but larger 
empty regions. We note that this configuration easily embeds isometrically 
in a square. 

Besides the discrete geometry background, we are also motivated by a 
problem in statistical physics, and in fact this was our original motivation. 
Hard discs in a box is one of the most basic mathematical models for mat- 
ter. Here one considers the space of all possible configurations of n non- 
overlapping discs of radius r in a unit square. This space has a natural 
probability measure, which we will not discuss here. But conceptually one 
could think of noninteracting particles, or particles where attractive forces 
are so much weaker than repulsive forces that the repulsive forces domi- 
nates, and so in the limit making a constraint to non overlapping particles 
makes sense. See the background and references in Section 4 of [jSl]. It is of 
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Figure 2. Another planar stable configuration. This also 
gives examples of stable configurations in a square, for a 
sequence of increasing n and r = 1^(1/ v^)- 

theoretical interest to try to understand basic properties of phase diagrams, 
even in this simple mathematical model. 

Countless hours of computer simulations have been run with hundreds of 
discs, and various densities of r. (For example, see Chapter 3 of [,11.1 ). It 
has been observed that there is a phase transition when nr'^n ^ 0.70. In 
particular if nr'^n <^ 0.70 then the discs seem to be behave more "liquid," 
and if nr'^n ^ 0.70, they are more "solid," according to various statistics 
measured, indicating action at a distance, and so on. To sample a point in 
Config(n; r) uniformly, the Metropolis algorithm and its variants are used. 

The most basic implementation of Metropolis is to start with a particular 
configuration of discs, then iterates the following algorithm: choose one of 
the n discs uniformly randomly, and then perturb it randomly, by moving 
it uniformly within a smaller disc of radius e. If the proposed new position 
of the disc does not go out of the box, and does not overlap with any other 
discs, then accept the move. Otherwise, reject it. Repeat this process as 
necessary. 

In order to make the "as necessary" precise, Diaconis, Lebeau, and Michel 
recently established useful bounds on the mixing time of the Metropolis al- 
gorithm in general, and then applied their results to the specific case of hard 
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discs in a box [6]. In order to meaningfully discuss mixing time questions, 
they must require that the discs be sufficiently small that one can move 
from every configuration to every other configuration, via small moves by 
one disc at a time (i.e. that the Metropolis Markov chain is irreducible), and 
they give an upper bound on the disc radius that guarantees that this is the 
case. Let Config(ri; r) denote the space of all possible configurations of n 
non- overlapping discs of radius r in a unit square. A stable configuration 
represents an isolated point in Config(n; r). 

Theorem 1.1 (Diaconis, Lebeau, and Michel). There exists a constant a 
such that the Metropolis Markov chain on Config{n; r) is irreducible when- 
ever r < a/n. In particular, there are no stable configurations in this 
range. 

One might wonder if this result can be extended to larger r. Is it still true, 
for example, if r = 0{^/n)l However, just as in Fejer Toth's boundary-free 
version of this problem, the answer is no. The result of Diaconis, Lebeau, 
and Michel is is essentially best possible. 

Theorem 1.2. There exists a constant (3 such that for sufficiently large n 
and r > (3/n, there exist stable configurations ofn discs of radius r. 

This is our main result. We prove Theorem 1 1.2 1 in the next section. 

2. Sparse stable configurations 

A positioning of n nonoverlapping discs of radius r in the unit square 
[0, 1]^ is a configuration. We say that a configuration is stable if each disc 
is held in place by its neighbors and the walls. Our main result is the fol- 
lowing. 

We first describe a construction of a one-way infinite "bridge" in the 
plane, due to Boroczky We follow closely Pach and Sharir's descrip- 
tion in [flill . For convenience of notation, all of the discs for now are unit 
radius; they can be rescaled as necessary later to fit in a box. First, the 
construction is symmetric about the a;-axis, so we only describe discs with 
their centers on our above the axis. Let ai = (0, 2 + ^3), hi = (0, ^3), 
Ci = (1, 0), as in Figure [51 Let f{x) be strictly convex function defined for 
X > 0, with /(O) = 2 + -\/3 and lima.^oo f{x) = 2^/3. Let C be the curve 
C = {{x, f{x)) I a; > 0}, and let ai, 02, eta, . . ., the unique sequence of 
distinct points on C such that d(ai, flj+i) = 2 for i = 1, 2, . . .. 

Now set 62 to be the unique point to the right of 02 such that d{b2, 02) = 2, 
and d{b2, ci) = 1. Then set 02 to be the point on the x-axis to the right of 62 
with d{c2,b2) = 1. Inductively define fej+i to be the point to the right of Ui 
such that d{hi+i, Oj+i) = 2 and d{hi^i, Ci) = 2. Similarly, set q+i to be the 
point on the x-axis to the right of such that d{ci^i, 6^+1) = 2. 
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Figure 3. Boroczky's bridge. (After Figure 9.2, p. 282, in UM-) 



The beginning of this construction is illustrated in the left side of Figure 
m The details that the construction results in a well-defined, stable config- 
uration (except along the four discs on the left) can be found in Boroczky's 
original paper [2] . He needed a second idea to get a stable configuration and 
disprove Fejer Toth's conjecture. He constructed a junction where three of 
these bridges could meet and hold all the loose discs on the end in place. 
Boroczky's junction is shown in Figure |4l 

We first use a slight modification of Boroczky's bridge construction, due 
to Pach and Sharir [,141 . For e > 0, replace f(x) above by the strictly convex 
function 

:= (l + e)/(x)-e/(0). 

We have /^(O) = /(O) = and lim^_^oo fe{x) < 2-\/3, so the sequence 

no longer continues indefinitely. By the intermediate value theorem, by 
varying e we can insure that for any arbitrarily large N, some is the last 
well-defined point, and its x-coordinate is one more than the x-coordinate 
of a^v. 

Then let / be the vertical line through b^, as in Figure [51 and complete 
the configuration using / as a line of symmetry, as illustrated. This gives an 
arbitrarily long symmetric bridge. 
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Figure 4. Boroczky's junction. The dark shaded discs are 
stable, and three copies of the bridge will overlap along the 
lightly shaded discs. This gives a zero-density stable config- 
uration in the plane. 

The second piece of our construction for the square is the comer "junc- 
tion" piece in Figure [6l The construction is easily verified to exist. If the 
center of the bottom left disc has coordinates (0, 0), then the remaining five 
discs in the bottom left quarter of the square have coordinates: 

(0, 2), (2, 0), (2 + V3, 1), (1, 2 + ^3), and (1 + ^3, 1 + Vs). 

Then by taking four copies each of the junction and bridge, and arranging 
them so that they overlap along their shaded outermost disks, in Figure Ul 
the resulting configuration is stable. 

3. Future directions 

The configurations constructed in this article are stable, meaning that no 
single disc can move. This is important to know, particularly from the point 
of view of probability and statistical physics. But is it possible for several 
discs to move simultaneously? Given that one can continuously deform 
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Figure 5. Pach and Sharir's symmetric version of 
Boroczky's bridge. 




Figure 6. The comer junction. 
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Figure 7. The main construction. The white bridges can 
be stretched arbitrarily long, so there are stable configura- 
tions in the square with arbitrarily large n and r = 0{l/n). 



the strictly convex function f{x) in the main construction in Section [2l and 
if one is careful preserve the combinatorial structure (tangency graph), it 
would seem that there is a small amount of flexibility and that one can 
slightly perturb many disks at once. How much freedom is there, though? 
In particular, is it possible to freely permute the discs? 

The set of all configurations of n hard discs in a unit box, Config(n; r), 
inherits a subspace topology from M^". Persi Diaconis points out in his sur- 
vey article [5] that, "very, very little is known about the topology of these 
spaces." The question above about moving multiple discs simultaneously is 
most naturally phrased in terms of path components of Config(n; r). Con- 
sidering Config(n; r) as a topological space allows us to ask many new 
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Figure 8. A stable configuration of 5 discs that illustrated 
non-monotonicity of the property "path-connectedness." 
This configuration is an isolated point in Config(ri; r), but 
for slightly larger radius, computer experiments suggest that 
Config(5; r) is path connected. 

questions, as well. For example, we could consider path components un- 
der various restrictions (Lipshitz, smooth, etc.). If r is sufficiently small, 
it would seem that Config(ri; r) is not only path connected, but homo- 
topy equivalent to the classical configuration space Config(n, M^), a space 
well-studied in algebraic geometry and algebraic topology What is 
the largest R = R{n) such that Config(n; r) is homotopy equivalent to 
Config(n,R2) forr < Rl 

One might guess that once r is sufficiently small that Config(?2; r) is path 
connected, that if r' < r then Config(r;,; r') is path connected as well. How- 
ever, this is not the case. For example, the configuration in Figure [8] is 
stable, and in fact, it can be shown, that it is a single point component in 
Config(n; r). But for slightly larger radius, it would seem that Config(n; r) 
is path connected. 

On the subject of non-monotonicity, we note the work of Brightwell, 
Haggstrom, and Winkler , where hardcore and Widom-Rowlinson models 
are studied on graphs. In [4] the property of having a unique Gibbs mea- 
sure, perhaps somewhat analogous to Config(n; r) being path connected, is 
shown to be non-monotone in the underlying parameter. 

The topological questions may be interesting for their own sake, but they 
may also have implications for understanding the statistical physics. For ex- 
ample, although we have showed that Metropolis is not irreducible on this 
problem, even for reasonably small discs, some might complain that the 
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Stable configurations only comprise a set of measure zero in Config(n; r). 
Suppose that one could show that Config(n; r) is disconnected. Then one 
can use a compactness argument to show that Config(n; r') is still discon- 
nected for some r' < r. In fact for some r' < r, 

|7ro(Config(n;r))| < |7ro(Config(?2; r'))| 

, where 7ro(X) denotes the set of path components of X. There is a natu- 
ral inclusion map i : Config(n; r) ^ Config(n; r'), and every element of 
i(7ro(Config(n, r)) has positive measure in Config(n; r'). 

For example, in Figure [8l one can slightly shrink the discs, and there are 
still at least 4 x 5! + 1 = 481 path components, each with small positive 
measure. 

Other discrete geometry questions naturally arise. Can an analogue of 
the construction in this article be carried out in higher dimensional cubes? 
What about for an arbitrary smoothly bounded convex body in the plane? 
(The result of Diaconis, et. al [6] is much more general than we stated 
earlier. In particular. Theorem 11.11 holds for general convex domains with 
reasonably nice boundaries.) 

The most important question in this area, from the point of view of statis- 
tical physics, is to give a mathematical explanation for the phase transition 
observed experimentally to occur when vrr^n ^ 0.70. To prove that any 
phase transition occurs at all in the hardcore model seems to be an open 
problem, although we note [3], where a a phase transition is proved to exist, 
in a similar hardcore model, but where the particles are nonconvex . 
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